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ABSTRACT 

We derive a semi-analytic solution for the structure of conduction-mediated 
transition layers above an X-ray illuminated accretion disk, and calculate 
explicitly the X-ray line radiation resulting from both resonance line scattering 
and radiative recombination in these layers. The vertical thermal structure of 
the illuminated disk is found to depend on the illuminating continuum: for a 
hard power law continuum, there are two stable phases connected by a single 
transition layer; while for a softer continuum, there may exist three stable 
phases connected by two separate transition layers, with an intermediate stable 
layer in between. We show that the structure can be written as a function of 
the electron scattering optical depth through these layers, which leads to unique 
predictions of the equivalent width of the resulting line radiation from both 
recombination cascades and resonance line scattering. We find that resonance 
line scattering plays an important role, especially for the case where there is no 
intermediate stable layer. 

Subject headings: accretion disks - galaxies: nuclei - radiative transfer - resonant 

scattering - X-rays:spectra 
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1. INTRODUCTION 

X-ray reflection off the surface of cold disks in active galactic nuclei (AGN) and galactic 
black holes (GBHs) has been an active field of research since the work of [Lightman fc White 



|(1988|) . In early studies, the illuminated material was assumed to be cold and non-ionized 



( [Lightman fc White 1988t [George fc Fabian 1991|) . It was soon realized, however, that 



photoionization of the disk can have a great impact on both the reflected continuum and 
the iron fluorescence lines. Detailed calculations were then carried out by [Ross fc I'abian 



(TO5[ ); [Matt, ^'abian fc Ross (1993[ ); [Zycki et al. (ig^ ); [Ross, ^'abian fc Brandt (19961) and 



Hoss, Fabian fc Young (1999 ). However, in all of these papers, the density of the illuminated 



material was assumed to be constant along the vertical direction. This assumption applies 
only to the simplest version of radiation-dominated Shakura-Sunyaev disks ([Shakura &j 



punyaev 1973|) , and only for the portion where viscous dissipation is the dominating heating 
process. For the surface layers, however, photoionization and Compton scattering are the 
major heating sources. Therefore the approximation of constant density is not appropriate. 
Moreover, thermal instability allows the coexistence of gas at different phases. These 
different phases have very different temperatures, and hence different densities to keep the 
gas in pressure balance. 



Recently INayakshin, Kazanas fc Kallman (2000 ) relaxed the simplifying assumption 



of constant gas density. They determined the gas density from hydrostatic balance solved 
simultaneously with ionization balance and radiative transfer. They made an important 
observation that the Thomson depth of the hot coronal layer can have great influence on 
the X-ray reprocessing produced by the deeper, and much cooler disk. In order to simplify 
the calculation of the vertical structure, though, they ignored thermal conduction and the 
effects of transition layers between the different stable phases. A discontinuous change in 
temperature was allowed whenever an unstable phase was encountered. They argued that 
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such transition layers are of little importance because their Thomson depths are negligibly 
small. However, without taking into account the role of thermal conduction, their method 
of connecting two different stable layers is rather ad hoc. Moreover, even though the 
Thomson depths of these transition layers are small, it does not guarantee that the X-ray 
emission and reflection from such layers are negligible. Because the temperature regime 
where the transition layers exist is not encountered in the stable phases, some of the most 
important lines can have appreciable emissivity only in these layers. Also, since resonance 
line scattering has much larger cross section than Thomson scattering, the optical depths 
in resonance lines can be significant. 

Including thermal conduction in the self-consistent solution of the vertical structure 
presents a serious numerical challenge. The difficulties are due to the coupling between 
hydrostatic balance, radiative transfer and heat conduction. [Zel'dovich fc PikeUner (1969|) 



first studied the phase equilibrium of a gas heated by cosmic rays and cooled by radiation. 
They found that taking into account heat conduction in the boundary layer allows one 
to obtain a unique solution of the stable equilibrium. [Rozahska (1999|) calculated the 
full temperature profile for a Compton- heated corona, and Rozahska fc Czerny (2000a ) 



calculated the static conditions of the plasma for different forms of heating and cooling. 
But they did not include much discussion of the spectroscopic signatures resulting from the 
derived vertical structure. 

In this paper, we first calculate the temperature structure in the layers above the 
accretion disk, then calculate the emission lines via radiative recombination (RR) and 
reflection due to resonance line scattering from the derived layers. Certain illuminating 
continua spectra allow more than two stable phases to coexist, with two transition layers 
connected by an intermediate stable layer. For the transition layer, since the Thomson 
depth is small, the ionizing continuum can be treated as constant; and since its geometric 
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thickness is smaller than the pressure scale height, the pressure can be treated as constant 
as well. We can thus obtain semi-analytic solution of the temperature profile by taking into 
account thermal conduction. For the intermediate stable layer, its thickness is determined 
by the condition of hydrostatic equilibrium. In our model, the normally incident continuum 
has a power-law spectrum with an energy index of a = 1. We also assume a plane-parallel 
geometry and that the resonance line scattering is isotropic. 

The structure of this paper is as follows: In §H we discuss the existence of the thermal 
instability and compute the thermal ionization structure of the transition layers; in §^ 
we calculate the recombination emission lines and the reflection due to resonance line 
scattering; in §^ we summarize the important points of the calculations, the validity of 
various approximations made in the calculations, and the detectability of the recombination 
emission and reflection lines. 



2. THERMAL INSTABILITY AND TRANSITION LAYERS 

2.1. The "S curves" 

The vertical structure of an X-ray illuminated disk at rest is governed by the equations 
of hydrostatic equilibrium and of energy conservation 



In the first equation, F is the force density due to gravity and radiation pressure. The 
dependence of the force on the plasma density is included explicitly through the hydrogen 
density Hh^x). In the second equation, a time independent state is assumed, /t(T) is 
the thermal conductivity, and ^{T, P, Fx{E)) is the net heating rate depending on the 
gas state and the incident flux Fx{E) (differential in energy). We neglect the effects of 





(2) 
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magnetic field and adopt the Spitzer conductivity appropriate for a fully ionized plasma, 
k(T) = 5.6 X 10-^T5/2 

erg cm^ s ^ K ^ (Spitzer 1962|) . We have used the classical heat flux, 



Qciass = — k(T)VT, in Equation (||) because the electron mean free path is short compared 
to the temperature height scale. 

Since the continuum flux may change along the vertical height, in principle, the above 
two equations must be supplemented by an equation for radiative transfer. A self-consistent 
solution of such equations is difficult to obtain. In the following, we invoke a few physically 
motivated approximations, which make the problem tractable. 

First, in thermally stable regions, the gas temperature is slowly varying, the heat 
conduction term in the energy balance equation can be neglected. Therefore, the 
temperature can be determined locally with the condition ^{T, P, Fx{E)) = 0. It is well 
known ( p:lrolik, McKee fc 'I'arter 1981| ) that the dependence of $ on the gas pressure P 
and the illuminating continuum Fx{E) can be expressed in the form of $ = —nlAnet{T,r.), 
where is the electron density, Anet(T, H) is the net cooling rate per unit volume, and H is 
an ionization parameter defined by 

where is the total fiux of the continuum, and c is the speed of light. In Figure |I], we 
show the local energy equilibrium curve T versus H, at Anet(T, S) = calculated with 
the photoionization code XSTAR ( |Kallman fc Krolik 1986| ). These curves are commonly 
referred to as "S curves" due to their appearance. The illuminating continuum is assumed 
to be a power law with energy index a = 1. The solid line labeled with "S-curve 1" 
corresponds to a low energy cutoff at 1 eV and a high energy cutoff at 150 keV, while the 
dashed line "S-curve 2" corresponds to a high energy cutoff at 200 keV. The choice of such 
incident spectra is based on their common appearance in many AGNs and GBHCs. The 
region is thermally "unstable" where the "S-curve" has a negative slope, and "stable" where 



the slope is positive, as indicated in Figure |l]. In the thermally unstable regions, we have 

< 0, (4) 



df 



Anet=0 



where the derivative is taken while the energy balance is satisfied, i.e., Anet = 0. This 
condition was shown ( |Krolik, McKee fc Tarter 1981| ) to be equivalent to the instability 
condition discovered by [Field (1965D 

/ dAnet'^ 



dT 



< 0. 



(5) 



EDITOR: PLACE FIGURE |l| HERE. 



2.2. The Temperature Profiles 

2.2.1. The Transition layers 

The thermal instability allows the gas to coexist at different phases. The gas 
temperature may change by orders of magnitude over a geometric thin region whenever an 
unstable phase separates two stable ones. This results in enormous temperature gradients 
and heat conduction. Therefore the heat conduction in the energy balance equation should 
be included in such transition layers between stable phases. On the other hand, the 
thicknesses of these transition layers are usually smaller than the pressure scale height, so 
one can safely treat the gas pressure as constant in these regions. Moreover, the continuum 
radiative transfer can be neglected because the Compton optical depth is found to be small. 
The vertical structure of the transition regions is then solely determined by the energy 
balance equation with heat conduction. The Thomson optical depth drg = n^aTdx of such 
regions is readily estimated by 
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where ctt is the Thomson scattering cross section. 

The transition layer solution is not arbitrary under the steady-state conditions, i.e., 
where there is no mass exchange between the two stable phases which the transition layer 
connects. A similar problem in the context of interstellar gas heated by cosmic rays was 
well studied by [Zerdovich fc PikeFner (1969| ). We follow their procedure here and define 
y = K,{T)dT/dTe, in order to rewrite Equation (0) in the form: 

dy^ 2/€(T)A„ei(T,S) 



dT 4 

A steady-state requires vanishing heat fiux at both boundaries of the transition layer, or 
y |r=Ti -y \t=T2= = -2 / 2 dT, 



(7) 



where Ti and T2 are the temperatures of the two stable phases which are connected by 
a transition layer. This condition determines a unique ionization parameter H for the 
transition layer in question, and the integration of Equation (|^) along the vertical height 
gives the detailed temperature profile as a function of optical depth. 

If the disk does not realize the steady-state solution, there are additional enthalpy 
terms in Equation which require that there be mass fiow through the transition region 
- i.e. the cool material in the disk evaporates, or the hot material in the corona condenses 
( IMeyer fc Meyer- Hofmeister 1994t ^ullemond 1999| ; [Rozahska fc Czerny 2000b|) . Physically, 
this corresponds to a movement of the transition layer up or down through the vertical disk 
structure. However, since the density increases monotonically toward the center of the disk, 
this "motion" stops where the transition layer reaches the steady state value of S. Thus, in 
the absence of disk winds or continuous condensation from a disk corona, the steady state 
solution should generally be obtained. 
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2.2.2. The Intermediate Stable Layer 

There is a complication for the "S-curves" shown in Figure 0, because in each curve 
there exist two unstable regions and therefore there should be two transition layers. For 
"S curve 1" , Condition ^ can be met for both transition layers with resulting ionization 
parameters Si < S2, where Si and S2 are associated with the transition layer which 
connects to the lowest temperature phase and highest temperature phase, respectively. For 
"S curve 2", the resulting ionization parameters for two transition layers, however, satisfy 
Hi > S2. Such a situation is unphysical, where the ionization parameter of the upper 
transition layer is smaller than that of the lower one, because in the context of accretion 
disks, the upper layer receives more ionizing flux and has lower pressure. So in practice, 
for "S curve 2" the intermediate stable region is skipped and a transition layer connects 
the lowest temperature phase to the highest temperature phase directly. The ionization 
parameter of this transition layer is determined the same way by applying Equation (H). 

There is then an intermediate stable layer, of nearly uniform temperature, in between 
these two transition layers for "S curve 1", as indicated with BC in Figure 0. The thickness 
of the intermediate layer should in principle be obtained by solving the coupled equations of 
hydrostatic equilibrium, energy balance, and radiative transfer. However, unlike the stable 
phase at the disk base or that of the corona, which may be Compton thick, this intermediate 
stable layer is generally optically thin, because its optical depth is restricted by the 
difference between Si and H2. Furthermore, the temperature in this layer is slowly varying, 
and therefore heat conduction can be neglected. We shall make another approximation 
that the variation of the force density F in the hydrostatic equation can also be neglected. 
This may not be a good approximation. However, since our main purpose is to investigate 
qualitatively the effects of an intermediate stable layer, such a simplifying procedure does 
capture the proper scaling relations of the problem, and has the advantage of less specific 
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model dependence. Writing Equation (|l]) in a dimensionless form and parameterizing the 
force factor by a dimensionless parameter A = cF/aTF^, we obtain 



1 (is , 



The parameter A defined here is identical to the gravity parameter in [Nayakshin, Kazanas 



fc Kallman (20001 ) in the absence of radiation pressure. The integration of this equation 



from Hi to H2 gives the Thomson depth of the intermediate stable layer. 

Assuming radiation pressure can be neglected, The force factor at a given radius R of 
the disk can be estimated as 

lHAnGMm„c 
^ = J-B J (10) 

where / is the luminosity of the continuum source in units of the Eddington limit, H is 
the half thickness of the disk at radius R, M is the mass of the central source, rrip is the 
proton mass, and Ledd is the Eddington luminosity. In this estimate, we have assumed 
that Fx = L/AtiR^. This gives A the value: 

A^1a\1. (11) 

For a typical thin disk, as those present in AGN and black hole binaries, one expects 
H/R ~ 0.01. If the luminosity is sub Eddington, / ~ 0.01 as in most AGNs, A is of order 
unity. However, since the disk surface may not be normal to the continuum radiation, 
Fj. may be only a fraction of the value assumed above, which increases A by one or two 
orders of magnitude. On the other hand, as the source approaches the Eddington limit, A 
may become smaller than 1. Therefore, we expect A to be in the range of 0.1 - 10. The 
exceptional cases of much smaller and larger A are discussed in §^ 

The temperature profiles and optical depths of the transition layers and possible 
intermediate stable layer are shown in Figures |^. Three labeled curves in solid lines 
correspond to "S-curve 1" (in Figure |I|) with different values of the force factors A = 10, 
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1, and 0.1, respectively. The dashed curve corresponds to "S-curve 2", where there is no 
intermediate stable layer. For each of the solid curves, three layers are clearly seen, with two 
transition layers being connected by an intermediate stable layer, as illustrated for the case 
of y4 = 1. The smaller A produces a more extended intermediate stable layer as expected. 



EDITOR: PLACE FIGURE | HERE. 



3. X-RAY EMISSION AND RESONANCE LINE SCATTERING 

3.1. X-ray Emission From Radiative Recombination (RR) 

Because the stable phase with the lowest temperature is almost neutral, and the stable 
phase at the highest temperature is almost fully ionized, they are not efficient in generating 
X-ray line emission, except for iron fluorescence lines from the neutral material. Only the 
transition layers and the intermediate stable layer are expected to emit discrete lines in 
the soft X-ray band. In a photoionized plasma, the temperature is too low for coUisional 
excitation to be an important line formation process. Instead, radiative recombination 
(RR) followed by cascades dominates the line emission. The flux of a particular line can be 
written as 

Fem = J n,niEiAi{T)dx, (12) 

where Ui is the density of the ion before recombination, Ei is the line energy, and A;(T) is 
the line emissivity defined as in [Sako et al. (1999| ) . 

In ionization equilibrium where the ionization rate is equal to the recombination rate, 
we have UeUiai = / ^""^^ ai^i[E)dE , where is the recombination coefficient of ion 

nj_i is the number density of ion i — 1, Fx{E) is again the monochromatic incident 
flux (differential in energy), and crj_i(-E) is the photoionization cross section of ion i — 1. 
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Defining the branching ratio pi{T) = Equation (|T2|) can be rewritten as: 

= EiJ ^^"-^^dE J MT)fU^, T)dT,, (13) 

where /j_i(T, S) = ^^^^ is the fractional abundance of ion i — 1 with respect to the electron 
density. 

As indicated, pi and depend only on temperature T and ionization parameter 
S, which are both functions of optical depth Tg. For convenience, we further define the 
"emission equivalent width" EWi^ = Y^f^- Then if h{E) = -f^^^, which depends only on 
the shape of the incident continuum, EWl' can be written as: 

= EiJ ^^^^dE J p,(r)/._i(T,H)rfre. (14) 

Therefore the emission equivalent width EWf is independent of the density of the 
medium and the incident flux. It is a unique function of the structure deduced in Section ^ 
All other variables depend only on T and S. The numerical values of pi{T) were provided 
by D. Liedahl (private communication), and were calculated using the models described 
in ^ako et al. (1999[) . The values of /i_i(T, S) were computed using XSTAR ( [Kallman fc 
Krolik 19861) . 

In Figure ^, we plot the spectra of the recombination emission within the 0.5 - 1.5 
keV band with a spectral resolution ~ 2 eV, which is close to the spectral resolution of 
the grating spectrometers on Chandra and XMM-Newton. The top panel corresponds to 
the case without an intermediate stable layer, and the bottom panels corresponds to the 
case with an intermediate layer for A = 10,1 and 0.1, respectively. For clarity, from top to 
bottom, the flux in each panel is multiplied by a factor indicated in each panel. It appears 
that the existence of an intermediate stable layer enhances the emission in this energy 
band. This is not surprising since the ions that are responsible for these lines have peak 
abundances at temperatures close to that of the intermediate stable layer. In all cases. 
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the equivalent widths (EWs) of the emission hnes are less than 1.0 eV, with respect to 
the ionizing continuum. The strongest lines are the hydrogen-like and helium-like lines of 
oxygen, with EWs approaching several tenths of an eV. Hydrogen-like and helium-like lines 
from iron outside the plotted energy band are somewhat stronger, with EWs reaching a few 
eV. We note that our low equivalent width values conflicts with those derived by [Nayakshin^ 
Kazanas fc Kallman (200"(]| ), who found some lines with equivalent widths as high as 30 eV. 
However, since they did not consider the appropriate locations for the transition regions in 
S-space, their intermediate stable layer subtended a much larger optical depth than we find 
here. Naturally with a thicker layer, they found larger equivalent widths. 

EDITOR: PLACE FIGURE | HERE. 

3.2. X-ray Reflection From Resonance Scattering (RS) 

Emission from RR is not the only line formation process in the transition layers and the 
intermediate stable layer. Due to very large cross sections in resonance line scattering, the 
reflected flux in these lines may be significant. With the computed thermal and ionization 
structure, the column density in each ion and the optical depth in all resonance lines can 
be calculated straightforwardly. 

The cross sections for resonance line scattering depend on the line broadening. We 
assume thermal Doppler effects as the only mechanism. Although the gas temperature is 
a function of depth, we calculate the line width for a temperature where the abundance of 
each ion peaks as an average, and assume that the resonance scattering cross sections are 
uniform along the vertical direction. In terms of absorption oscillator strength /, this cross 
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section of the resonance line scattering a^is can be written as 

<yRS = ^/xT' 

where me is the electron mass, e is the electron charge, A is the wavelength of the line and 
AA is the average line width in wavelength. Under the assumption of thermal Doppler 
broadening, 

AA = 2.35A./^, (16) 

where Tm is the temperature at which the ion abundance peaks, and rrii is the ion mass. 
The resonance scattering optical depth trs for a line from ion i can be estimated as 

trs = (TRsNi, (17) 
where Ni is the column density of the ion. 



The radiative transfer in the line is a complicated issue ( Pumont, Abrassart fc Collin 



2000|) . A full treatment is beyond the scope of this work. However, since we are only 



interested in a reasonable estimate of the reflected line flux, a simple approach may be 
adopted. We assume the resonance line scattering is isotropic and conservative and neglect 
the polarization dependence. Under such conditions, the reflection and transmission 
contributions by a plane-parallel slab of flnite optical depth r have been solved by 
Uhandrasekhar (196CI|) . For normal incident flux F^, the angle dependent reflectivity is 



where /i = cos^, r(/i) is the reflectivity at /i, / is the reflected intensity, is the incident 
flux, and S{t, fj,) is the scattering function deflned as 

(1 + -)S(r, /i) = X(/i)X(l) - F(/i)r (1), (19) 
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where X(/i) and Y{^) are two functions that satisfy the following integral equations: 

= 1 + ^ ^'^^[X{^^)Xi^^')-Yi^^)Yi^,')W 
2 Jo /i + /i' 

Yi^^) = 6-^/^^ + ^ ['^-^[Yif^)Xif^')-Xif^)Yifi'W. (20) 
2 Jo /i — /i 

The solutions of these equations may be obtained by an iterative method with the starting 
point X{ii) = 1 and Y{ix) = e~'^^^. The angle integrated reflectivity Rt can be calculated as 

Rt = 27r [\{fx)fxdfx, (21) 
Jo 

and is shown in Figure | as a function of the resonance scattering optical depth tjis- 



EDITOR: PLACE FIGURE | HERE. 



The reflected flux in a line can be written as 

Fr = RtF,{E)AE, (22) 

where AE is the line width in energy. Similarly to the "emission equivalent width" EWf, 
we define a "reflection equivalent width" EWf = "p"^^? which results in: 

EWl = (23) 

This "reflection equivalent width" from our numerical results is a few tenths of an eV 
for strong resonance lines, similar to that of the recombination emission lines. 

In Figure |], we plot the spectra of the resonantly scattered lines in the energy band 
0.5 - 1.5 keV with a spectral resolution ~ 2 eV. The top panel corresponds to the case 
without an intermediate stable layer, and the bottom panels corresponds to the case with 
an intermediate layer for A = 10,1 and 0.1, respectively. For clarity, from top to bottom, 
the flux in each panel is multiplied by a factor of 100. We see that the equivalent widths of 
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the reflected lines are notably enhanced when there is an intermediate stable layer, but not 
as significantly as for the recombination emission lines. This is because the optical depths 
of many strong lines become much larger than unity, and the refiection is saturated. If there 
are broadening mechanisms other than thermal Doppler effects, such as turbulent velocity, 
the reflected intensity can be further enhanced. 

EDITOR: PLACE FIGURE | HERE. 



3.3. X-ray Spectrum with RR Emission and RS Refiection 

In order to gain a crude idea of the relative importance of recombination emission and 
reflection from the transition layers and intermediate stable layer, we compare them to the 
"hump" produced by Compton scattering off a cold surface ( [Lightman fc White 1988|) . We 
use the Greens function obtained by [Lightman fc White (19881) to calculate the Compton 
reflection. This method was verifled to be accurate with a Monte Carlo procedure by [George 
fc Fabian (1991|) . In Figure we show the combined spectra including recombination 
emission and reflection lines from resonance line scattering, and the Compton reflection 
"hump" . 



EDITOR: PLACE FIGURE | HERE. 



4. SUMMARY 

We now summarize the most important conclusions that can be drawn from the 
calculations presented in this paper. We also discuss the detectability of the predicted line 
features. 
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1. The unique ionization parameters that characterize the steady-state solutions of 
the transition layers depend on the shape of the "S-curve". We have shown that 
two power-law illuminating spectra with different high energy cutoffs produce very 
different temperature profiles. The harder spectrum only allows one transition layer 
even though there are two unstable branches in the "S-curve", while the softer one 
allows two separate transition layers connected by an intermediate stable layer. This 
is due to the fact that the ionization parameter of the upper transition layer must be 
larger than that of the lower one, if they are to exist separately in a disk environment. 
The harder spectrum produces a turnover point of the upper branch of the "S-curve" 
at smaller H. Therefore the transition layer due to the upper unstable region joins 
the lower one smoothly without allowing the intermediate stable region to form. The 
turnover of the upper "S-curve" represents the point where Compton heating starts 
to overwhelm bremsstrahlung. The ionization parameter at which this point occurs is 
related to the Compton temperature of the continuum, S^p oc Tj^^"^ ( [Krolik, McKee 



& 'iarter 1981; [Nayakshin 2(J(J(J ). A harder spectrum has larger T/c", therefore the 



intermediate stable layer tends to disappear for hard incident spectra. 

2. Although the Thomson depths of the transition layers and possible intermediate 
stable layer are generally negligible, the X-ray emission lines from them may comprise 
the main observable line features, because the temperatures of these layers are 
inaccessible to the stable phases, and thus some of the important lines can have 
appreciable emissivity only in these layers. Due to the much larger cross sections 
for resonance line scattering, reflection due to resonance lines off such transition 
layers is also important. The strengths of reflected lines are at least comparable with 
those of the recombination emission lines when there is no intermediate stable layer. 
Because the appearance of the reflected line spectrum is different from that of the 
recombination emission spectrum, high resolution spectroscopic observations should 
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be able to distinguish these mechanisms. 

3. The justification of the assumption that the ionizing continuum does not scatter in 
the intermediate layer depends on the magnitude of the parameter A. The Thomson 
depth of this layer Xe is given by: 

re = ^{^ - (24) 

For the power-law continuum with high energy cutoff at 150 keV ("S curve 1"), 
Si = 2.82 and S2 = 3.14 from our numerical results. Therefore, = 0.036/^4. t^. 
is much less than unity as long as A is greater than 0.1. For smaller A, however, 
another effect comes into play. [Nayakshin, Kazanas fc Kallman (2000|) showed that 



the Thomson depth of the coronal layer (the stable phase with highest temperature) 
exceeds unity when the gravity parameter (identical to A defined here when the 
radiation pressure is neglected) is ~ 0.01. Therefore not much ionizing flux can 
penetrate this layer, and the reprocessing in the deeper and cooler layers can be 
neglected completely. As A becomes much larger than 10, the thickness of the 
intermediate stable layer is negligible compared to the transition layers. Therefore, its 
presence may be ignored. 

4. Since the recombination rate must equal the photoionization rate in the irradiated gas, 
recombination radiation is also a form of reflection - i.e. the line equivalent widths 
are independent of the incident flux. They depend only on the structure (T(r); H(r)) 
deduced from the hydrostatic and energy balance equations. 

5. The detectability of these recombination emission and reflection lines depends on 
whether the primary continuum is viewed directly. When the ionizing continuum is in 
direct view, our results show that the EWs of the strongest lines in the 0.5 - 1.5 keV 
band are at most a few tenths of an eV, slightly larger when there is an intermediate 
stable layer. 



The signal to noise ratio (SNR) in such a hne can be written as 



EW 



itEiF.y/' (a 



E 



) 



1/2 



SNR^ 




(25) 



where t is the integration time of the observation, Ei is the energy of the hne, Fj is 
the photon flux in the continuum, Agff and E/5E are the effective area and resolving 
power of the instrument, respectively. For a line at ~ 1 keV, with EW = leV, and 
with HETGS on board Chandra, we have SNR = 0.25./tF\. A typical Seyfert 1 



galaxy has a flux of 10"^'^ erg cm~^ s~^ in the energy band of 2-10 keV. Assuming a 
power law with energy index of 1, the photon flux at 1 keV would be ~ 0.05 cm~^ 
s~^ keV~^. For a reasonable integration time of 10 ks, we have SNR ~ 5. When the 
primary continuum is obscured as in Seyfert 2 galaxies, the EWs of the emission 
and reflection lines can be orders of magnitude larger, because the continuum at this 
energy region is absorbed severely, and the SNR can be greatly enhanced, making 
these lines observable. 
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Fig. 1. — The "S-curves" produced by two different incident ionizing spectra. The vertical 
hne indicates the unique solution of S which satisfies Condition two solutions for "S-curve 
1", Si = 2.82 and S2 = 3.14; and only one for "S-curve 2", S = 2.22. 
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Fig. 2. — The temperature profiles of the transition layers and the intermediate stable layer 
versus Thomson optical depth Tg. The sohd curves correspond to "S-curve 1" with different 
force factors A, the dashed fine corresponds to the "S-curve 2" . 
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Fig. 3. — The spectra of the emission hnes via RR in the transition layers and the 
intermediate stable layer. The spectral resolution is AE ~ 2 eV. The top panel corresponds 
to "S-curve 2" , while bottom panels correspond to "S-curve 1" with different parameter A. 
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Fig. 4. — The angle integrated reflectivity of a plane parallel slab with normal incident flux 
as a function of the resonance scattering optical depth trs- 
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Fig. 5. — The spectra of the reflection hnes due to resonance scattering in the transition 
layers and the intermediate stable layer. The spectral resolution is ~ 2 eV. The top panel 
corresponds to "S-curve 2", while bottom panels correspond to "S-curve 1" with different 
parameter A. 
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Fig. 6. — The combined spectrum. Green - emission lines via recombination, red - reflection 
lines due to resonance line scattering, and black - the Compton reflection hump. The 
spectral resolution is ~ 2 eV. The top panel corresponds to "S-curve 2", while bottom 
panels correspond to "S-curve 1" with different parameter A. 



